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We  explore  two-dimensional  self-assembly  of  tobacco  mosaic  viruses  (TMVs)  on  a  substrate-supported, 
fluid  lipid  monolayer  by  manipulating  the  electrostatic  interactions,  with  specific  focus  on  the  effects  of 
the  cationic  lipid  concentration  in  the  monolayer  and  the  presence  of  Ca2+  ions  in  the  surrounding  bulk 
solution.  The  TMV  assemblies  were  characterized  by  grazing-incidence  X-ray  scattering  and  atomic  force 
microscopy,  and  the  inter-particle  interaction  quantified  through  X-ray  scattering  data  analysis.  In  the 
absence  of  Ca2+  ions,  we  found  that  higher  charge  densities  on  the  lipid  monolayer  led  to  poorer  in-plane 
order,  which  may  be  attributed  to  faster  adsorption  kinetics,  due  to  the  surface  potential  that  increases 
with  charge  density.  At  the  same  time,  higher  lipid-charge  densities  also  resulted  in  weaker  repulsion 
between  TMVs,  due  to  partial  screening  of  Coulomb  repulsion  by  mobile  cationic  lipids  in  the  monolayer. 
The  lipid-charge  dependence  was  diminished  with  increasing  concentration  of  Ca2+  ions,  which  also  led 
to  tighter  packing  of  TMVs.  The  results  indicate  that  Ca2+  ions  strengthen  the  screening  of  Coulomb  repul¬ 
sion  between  TMVs  and  consequently  enhance  the  role  of  attractive  forces.  Control  experiments  involv¬ 
ing  Na+  ions  suggest  that  the  attractive  inter-TMV  interaction  has  contributions  from  both  the  van  der 
Waals  force  and  the  counter-ion-induced  attraction  that  depends  on  ion  valence. 

©  201 1  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

Self-assembly  into  ordered  structures  has  been  widely  studied 
because  of  the  interests  in  utilizing  controlled  self-assembly  for 
practical  applications  [1,2].  Self-assembly  processes  in  general 
are  complex  and  can  be  influenced  by  factors  such  as  the  proper¬ 
ties  of  the  building  blocks  themselves  (e.g.  monodispersity,  shape), 
the  interaction  between  them,  and  kinetic  properties  of  the  system 
(e.g.  mobility,  reversibility).  In  simplified  model  systems,  it  is  often 
possible  to  control  these  factors  individually  and  quantitatively 
characterize  the  resulted  structure  so  that  the  influence  of  the  con¬ 
trol  factor  can  be  understood.  For  instance,  the  interaction  between 
the  building  blocks  can  be  modified  by  increasing  the  screening  of 
the  electrostatic  repulsion,  by  introducing  attraction  through  cap¬ 
illary  [3-5]  and  depletion  [6,7]  forces,  or  by  surface-functionalizing 
the  building  blocks  so  that  they  would  selectively  interact  with 
each  other  [8,9]. 
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Recently  we  reported  on  two-dimensional  assemblies  of  tobac¬ 
co  mosaic  viruses  (TMVs)  and  the  use  of  X-ray  scattering  to  char¬ 
acterize  the  TMV-TMV  interaction  within  the  assemblies  [10].  In 
this  model  system,  the  negatively  charged  TMVs  are  adsorbed  on 
a  substrate-supported,  positively  charged  lipid  monolayer  at  a  bulk 
buffer  solution-substrate  interface.  The  electrostatic  interaction 
between  these  charged  rods  can  be  controlled  by  adjusting  the 
solution  properties  (e.g.  pH  and  ionic  strength),  while  the  fluidity 
of  the  lipid  monolayer  gives  the  particle  mobility  that  is  necessary 
to  assemble  into  ordered  structures.  Due  to  the  rod-like  shape 
(300  nm  long  and  18  nm  in  diameter)  of  TMV  particles,  the  assem¬ 
blies  consist  of  TMVs  that  are  oriented  along  the  interface  and 
aligned  parallel  to  their  neighbors  (Fig.  1 ).  Therefore  we  can  model 
these  assemblies  as  one-dimensional  fluids  that  are  characterized 
only  by  the  lateral  distances  between  neighboring  TMVs,  defined 
as  d,  while  neglecting  their  orientations.  An  attractive  attribute  of 
this  model  system  is  that  there  exist  theories  that  analytically  de¬ 
scribe  the  behavior  of  ID  fluids  for  several  types  of  inter-particle 
interactions  [11-13].  The  approximation  as  ID  fluid  therefore  al¬ 
lowed  us  to  fit  the  X-ray  scattering  data  from  the  TMV  assemblies 
based  on  these  theoretical  predictions.  This  approach  provided  a 
means  to  characterize  the  nature  of  the  inter-particle  interaction 
and  to  quantify  the  effects  of  external  controls  on  the  inter-particle 
interaction  and  ultimately  the  self-assembly  process. 
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Fig.  1.  (A)  Schematic  of  a  2D  TMV  assembly  and  the  sample  geometry  for  GISAXS  measurements.  The  qz  component  of  the  scattering  vector  is  along  the  sample  normal  while 
qr  is  parallel  to  the  sample  surface.  (B)  Representative  in-plane  scattering  intensity  profiles  along  the  surface  enhance  ridge  (blue  horizontal  dashed  line  in  A),  as  measured  at  a 
series  of  times  after  introducing  TMVs  into  the  bulk  solution.  The  in-plane  scattering  data  are  analyzed  by  model  fitting  (red  lines  in  B),  as  described  in  the  text,  to  characterize 
the  interaction  between  TMVs.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


The  methods  to  tune  the  electrostatic  interaction  between 
charged  building  blocks  include  varying  solution  pH  and  ionic 
strength,  which  respectively  control  the  net  charge  carried  by 
each  building  block  [14]  and  the  screening  of  Coulomb  repulsion 
between  them  [15].  The  electrostatic  interaction  between  TMVs 
can  also  depend  on  other  parameters.  For  instance,  our  previous 
study  provided  evidence  that  the  presence  of  Ca2+  ions  gave  rise 
to  attraction  between  TMVs,  based  on  the  observation  of  in¬ 
creased  TMV  packing  density  and  enhanced  structural  order  in 
the  assemblies.  Others  have  also  reported  similar  results  in 
which  highly  ordered  2D  structures  were  induced  by  multivalent 
ions  [16,17]. 

Here,  we  explore  this  model  system  further  by  examining  the 
effects  of  the  lipid  monolayer  charge  density  and  the  Ca2+  concen¬ 
trations  on  the  formation  of  ordered  TMV  assemblies.  Intuitively, 
one  may  expect  higher  monolayer  charge  density  to  promote 
TMV  adsorption  to  achieve  higher  TMV  density  and  perhaps  bet¬ 
ter  structural  order.  On  the  contrary,  the  present  measurements 
in  the  absence  of  Ca2+  show  that  the  structural  order  is  actually 
poorer  at  higher  lipid-charge  density,  likely  due  to  faster  TMV 
adsorption.  Moreover,  quantitative  analysis  of  the  grazing-inci¬ 
dence  small-angle  X-ray  scattering  (GISAXS)  data  demonstrates 
that  increasing  the  lipid-charge  density  reduces  the  range  of  elec¬ 
trostatic  repulsion  between  the  adsorbed  TMVs,  possibly  due  to 
partial  screening  by  the  mobile  cationic  lipids  (DOTAP)  in  the 
monolayer.  On  the  other  hand,  under  the  100  mM  Ca2+  solution, 
the  lipid-charge  density  is  found  to  have  little  effect  on  the 
assembled  TMV  structures.  Time-dependent  measurements  fur¬ 
ther  show  that  the  presence  of  Ca2+  ions  reduces  the  rate  of 
adsorption  of  TMV  onto  the  lipid  monolayer,  which  may  also  con¬ 
tribute  to  the  better  structural  order  in  the  TMV  assemblies.  The 
role  of  Ca2+  ions  is  further  clarified  by  control  experiments  using 
solutions  that  have  the  same  ionic  strength  but  contain  Na+  ions 
instead.  The  TMV  packing  density  under  the  Na+  solution  is  found 
to  be  intermediate  between  the  low-salt  (pH  buffer  only)  and 
Ca2+-ion  cases.  The  difference  between  the  low-salt  and  high-salt 
(Na+  or  Ca2+  ions)  results  indicates  that  the  van  der  Waals  attrac¬ 
tion,  which  does  not  depend  on  ion  valence,  likely  contributes  to 
facilitating  the  formation  of  better-ordered  structures  at  high  salt 
concentrations  where  the  screening  of  electrostatic  repulsion  is 
strong.  On  the  other  hand,  the  difference  between  the  Ca2+  and 
Na+  results  is  evidence  that  counter-ion  induced  attraction,  which 
depends  on  ion  valence,  also  plays  a  role  in  the  ordering  of  TMV 
assemblies.  Our  data  therefore  suggest  a  dual  role  of  Ca2+  ions  to 
provide  screening  to  the  Coulomb  repulsion  as  well  as  to  mediate 
attraction  between  TMVs. 


2.  Experimental  methods 

Sample  preparation  and  data  collection  were  performed 
following  previously  reported  procedures  [10,18],  as  briefly 
described  below.  The  lipid  monolayers  of  mixed  diol- 
eoylphosphatidylcholine  (DOPC)  and  l,2-dioleoyl-3-trimethylam- 
monium-propane  (DOTAP)  were  formed  on  OTS-coated  silicon 
substrates  in  a  Teflon  cell  by  incubating  the  substrate  for 
~30  min  under  a  solution  of  lipid  vesicles  prepared  by  sonication. 
After  incubation,  the  sample  cell  was  rinsed  thoroughly  using  the 
buffer  solution  to  ensure  that  free  vesicles  were  removed  and  that 
the  bulk  solution  had  the  desired  pH  and  salt  concentration. 
TMVs  were  then  added  into  the  sample  cell.  GISAXS  measure¬ 
ments  (schematically  shown  in  Fig.  1A)  were  performed  at  beam¬ 
lines  X21  and  X9  at  NSLS,  as  the  TMV  assemblies  formed  on  the 
lipid  monolayer.  The  structure  became  stable  after  6-14  h,  and 
the  sample  cell  was  rinsed  again  before  a  final  GISAXS  pattern 
was  recorded.  Atomic  force  microscopy  measurements  were  per¬ 
formed  on  an  Agilent  5500  AFM  in  the  non-contact  tapping  mode 
using  Olympus  TR400PSA  cantilevers  (nominal  stiffness  ~0.08  N/ 
m,  resonance  frequency  ~151<Hz  in  water).  The  samples  were 
prepared  following  the  same  procedure  but  in  the  AFM  liquid 
sample  cell.  Only  the  final  structure  was  measured  using  AFM, 
after  rinsing  the  sample  cell  to  remove  non-adsorbed  TMVs. 

While  the  interpretation  of  AFM  images  is  usually  straightfor¬ 
ward,  the  X-ray  scattering  data  require  non-trivial  analysis. 
Following  our  previous  study  [10],  we  limit  our  analysis  to  the 
in-plane  scattering  profile  near  qz  =  0,  and  approximate  the  TMVs 
as  parallel  rods  interacting  via  either  exclusively  hard-wall  repul¬ 
sion  (hard  rod  model)  or  hard-wall  repulsion  with  a  short-ranged 
attraction  (sticky  rod  model).  As  we  have  previously  described  in 
detail  [10],  this  analysis  approach  permits  structural  characteriza¬ 
tion  of  TMV  assemblies  in  terms  of  the  packing  fraction  r\  of  the 
rods,  the  hard-wall  radius  a  and,  in  the  case  of  sticky  rod  model, 
the  stickiness  parameter  a,  which  reflects  the  strength  of  the 
attractive  potential  well.  These  parameters  provide  the  basis  for 
evaluating  how  electrostatic  interactions  affect  structural  order 
within  the  TMV  assemblies.  Clearly,  the  hard  rod  and  sticky  rod 
interactions  are  only  approximations  of  the  overall  interaction  be¬ 
tween  the  TMVs,  which  may  contain  contributions  of  various  ori¬ 
gins  (e.g.  screened  Coulomb  repulsion,  van  der  Waal  force  and 
counter-ion-mediated  attraction),  as  we  will  discuss  below.  There¬ 
fore  the  dependence  of  these  simplified  parameters  provides  a 
semi-quantitative  description  of  how  the  TMV-TMV  interactions 
are  controlled  by  external  parameters.  But  the  numerical  value  of 
these  parameters,  for  instance  the  hard-wall  radius  cr,  may  not 
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directly  correspond  to  any  physical  characteristic  of  the  TMVs  (i.e. 
the  diameter  of  the  TMV). 

3.  Results  and  discussion 

In  general  structural  order  in  2D  depends  on  the  particle  density 
and  the  inter-particle  interaction,  both  of  which  can  in  principle  be 
influenced  by  the  electrostatic  interactions  between  neighboring 
TMVs  and  those  between  the  TMV  and  the  lipid  monolayer.  In  Sec¬ 
tion  1,  we  delineate  how  the  charge  density  in  the  lipid  monolayer 
affects  the  TMV-monolayer  and  inter-TMV  interactions  and  the 
assembly  of  TMVs  under  low-salt  conditions  (10  mM  MES  buffer 
only  to  maintain  the  pH  in  the  solution).  In  Section  2,  we  describe 
how  the  presence  of  100  mM  Ca2+  ions  alters  the  TMV-monolayer 
interaction  and  diminishes  the  assembly’s  dependence  on  the  cat¬ 
ionic  lipid  concentration.  In  Section  3,  we  examine  the  role  of  Ca2+ 
ions  further  by  discussing  the  results  of  additional  measurements 
at  a  series  of  Ca2+  concentrations  and  those  in  the  presence  of 
Na+  ions. 

3.1.  Effects  of  lipid  monolayer  charge  density  on  TMV  assembly  in  the 
absence  of  Ca2+ 

We  collected  GISXS  data  from  TMV  assemblies  under  a  10  mM 
MES  solution  at  pH  6.0  at  a  series  of  mole  fraction,  XTAP,  for  the  cat¬ 
ionic  DOTAP  lipid  in  the  monolayer.  The  measured  in-plane  scat¬ 
tering  intensity,  recorded  when  the  TMV  assemblies  have 
become  stable  at  6-14  h  after  the  injection  of  TMVs,  are  summa¬ 
rized  in  Fig.  2A.  Clearly,  the  details  of  the  qr-dependent  intensity 
modulations  in  the  data  vary  with  the  lipid  composition.  In  partic¬ 
ular,  as  the  DOTAP  fraction  is  reduced  (XTAP  from  1  to  0.048),  the 
broad  shoulder-like  feature  at  qr  <  0.03  A-1  shifts  to  lower  qr  and 
develops  into  a  better  defined  peak.  This  behavior  suggests  that 
the  decrease  in  the  lipid-charge  density  tends  to  both  reduce  the 
local  packing  density  and  extend  the  range  of  the  positional  corre¬ 
lations  between  aligned  TMVs.  In  the  case  of  the  neutral  lipid 
monolayer  (XTAP  =  0),  the  GISAXS  data  (not  shown)  displays  no 
intensity  modulations  that  are  indicative  of  inter-TMV  correla¬ 
tions,  suggesting  no  significant  TMV  adsorption  in  the  absence  of 
cationic  lipids. 


In  the  absence  of  multivalent  cations,  the  inter-TMV  interac¬ 
tions  are  likely  dominated  by  the  screened  Coulomb  repulsion. 
We  have  previously  shown  [10],  and  confirm  further  in  Fig.  2 A,  that 
the  model  fitting  (solid  lines)  based  on  hard-wall  repulsion  ade¬ 
quately  capture  the  salient  features  of  the  observed  in-plane  scat¬ 
tering  profiles  (symbols).  The  corresponding  best-fit  curves  for 
S(qxa )  and  the  pair  correlation  functions  g(x/cr)  are  shown  in 
Fig.  SI  in  Supporting  information.  Note  that  the  packing  fraction 
r]  is  based  on  the  effective  diameter  of  the  hard  rod,  which  is  great¬ 
er  than  the  TMV  diameter  because  we  approximated  the  longer- 
ranged  electrostatic  repulsion  between  TMVs  as  simple  hard-wall 
repulsion.  Therefore,  instead  of  q,  we  utilize  the  ratio  q/a 
(Fig.  2C)  as  a  better  measure  of  the  TMV  packing  density  within 
aligned  domains  (~l/d). 

3AA.  Cationic  lipids  contribute  to  the  screening  of  repulsion  between 
TMVs 

The  most  notable  trend  in  Fig.  2A  and  C  is  that  the  range  of  repul¬ 
sion  between  TMVs  decreases  as  the  cationic  lipid  fraction  in  the 
monolayer  is  increased.  Specifically,  as  the  DOTAP  fraction  is  raised 
from  5%  to  1 00%,  the  GISAXS-extracted  o  is  found  to  decrease  mono- 
tonically  by  a  factor  of  ~2  and  approach  the  TMV  diameter  of  1 8  nm 
(Fig.  2C).  This  observed  o  behavior  may  be  attributed  to  stronger 
screening  of  the  electrostatic  inter-TMV  repulsion  at  higher  lipid- 
charge  densities  due  to  the  DOTAP  molecules  acting  as  mobile  cat¬ 
ions  confined  in  the  plane  of  the  monolayer. 

The  monolayer  charge  density  is  of  the  order  of  e/A/  x  XTAP, 
where  A/  ~  0.7  nm2/lipid  [19],  or  0.07-1.4  e/nm2  for  the  DOTAP 
fraction  of  XTAP  =  0.048-1.  On  the  other  hand,  the  surface  charge 
density  on  TMV  (with  diameter  D  =  18nm  and  length 
L  =  300  nm)  is  estimated  to  be  QJ(nDL)  <  ~0.1  e/nm2,  based  on  pre¬ 
viously  reported  values  of  the  net  charge  per  TMV  in  aqueous  solu¬ 
tions  (Q  <  ~1800  e)  [10,20,21].  Therefore,  over  most  of  the  probed 
range  of  lipid  compositions,  the  charge  density  on  the  lipid  mono- 
layer  should  far  exceed  that  on  the  TMV  surface.  Consequently,  as 
we  will  discuss  further  below,  a  significant  amount  of  free,  mobile 
charges,  i.e.,  DOTAP  and  condensed  counter-ions,  are  expected  to 
be  present  at  the  lipid  monolayer-solution  interface,  which  should 
re-distribute  upon  the  adsorption  of  TMVs  to  partially  screen  the 
electrostatic  repulsion  between  the  TMVs. 
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Fig.  2.  The  in-plane  scattering  profiles  along  the  surface-enhanced  ridge  (see  Fig.  1A)  in  the  final  scattering  patterns  from  TMV  assemblies  adsorbed  to  substrate-supported 
DOPC/DOTAP  monolayers  under  (A)  Ca2+-free  MES  buffers  and  (B)  MES  buffers  that  contain  100  mM  Ca2+.  The  mole  fractions  of  DOTAP  (XTap)  in  the  corresponding  lipid 
monolayers  are  listed  beside  the  data.  The  data  are  shown  as  symbols,  shifted  vertically  for  clarity.  The  solid  lines  are  the  best-fits,  based  on  the  hard  rod  model  for  (A)  and  the 
sticky  rod  model  for  (B).  The  corresponding  radius  of  hard-wall  potential  o  and  packing  density  of  rods  r\\o  are  shown  in  (C). 
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This  type  of  screening  that  is  confined  to  the  interface  has  been 
examined  theoretically  for  the  simple  case  of  two  infinitely  long 
charged  rods  near  the  plane  of  a  2D  salt  solution  by  Menes  et  al. 
[22].  They  found  the  potential  energy  for  the  inter-rod  repulsion 
to  be  proportional  to  (2/d)2,  where  d  is  the  spacing  between  rods 
and  the  effective  screening  length  2  is  inversely  proportional  to 
the  charge  density  in  the  2D  salt.  This  is  qualitatively  consistent 
with  the  observed  trend  of  decreasing  TMV  spacing  with  increasing 
lipid-charge  density.  We  have  previously  estimated  the  net  charge 
of  each  TMV  to  be  ~-90e  [10].  This  result  was  obtained  by  the 
Monte  Carlo  simulation  based  on  the  screened  Coulomb  repulsion 
between  TMVs  where  the  screening  was  assumed  to  arise  from  free 
ions  in  bulk  solution  only.  The  additional  screening  by  mobile  ions 
in  2D  (i.e.,  DOTAP)  may  explain  why  the  above  estimated  TMV 
charge  is  much  lower  than  the  often-cited  value  of  ~-1800  e  per 
TMV  [20,21]. 

3 A. 2.  Density  of  adsorbed  TMV  weakly  depends  on  lipid  charge 

The  data  in  Fig.  2A  and  C  also  show  that  as  the  lipid-charge  den¬ 
sity  is  increased,  the  inter-TMV  packing  density,  as  measured  by 
the  ratio  r\\c,  only  increases  slightly.  These  GISAXS  results  are  cor¬ 
roborated  by  AFM  measurements.  Representative  AFM  images  of 
equilibrated  TMV  assemblies  for  XTAP  =  0.16,  0.4  and  1.0  are  shown 
in  Fig.  3,  together  with  the  corresponding  GISAXS  patterns.  The 
tendency  of  the  lipid-bound  TMVs  to  pack  closely  and  align  locally 
is  evident.  Determined  from  measurements  at  10  distinct  locations 
in  each  AFM  image,  the  average  inter-TMV  spacing  is 
39.5  ±  3.0  nm,  34.8  ±  2.2  nm  and  32.7  ±1.7  nm  for  XTAP  =  0.16,  0.4, 
and  1.0,  respectively.  The  ~20%  decrease  in  the  AFM-derived  in¬ 
ter-TMV  spacing  with  increasing  lipid-charge  density  is  consistent 
with  the  ~20%  increase  in  the  best-fit  values  of  the  GISAXS-ex- 
tracted  r\\o  over  the  same  lipid  composition  range  (black  symbols, 
Fig.  2C). 

The  observed  increase  in  the  packing  density  of  TMVs  with 
increasing  XTAP  agrees  with  the  expectation  that  the  adsorption 
of  the  anionic  TMVs  onto  a  cationic  lipid  monolayer  is  driven 
entropically  by  the  counter-ion  release  mechanism  [23-27].  The 
higher  the  lipid-charge  density  is,  the  greater  number  of  counter¬ 
ions  should  be  initially  bound  to  the  interface  and  hence  can  be  re¬ 
leased  into  the  aqueous  solution  through  the  adsorption  of  macro¬ 
ions.  Thus,  the  adsorption  of  TMVs  is  expected  to  increase  with  the 
lipid-charge  density.  This  is  consistent  with  the  observations  made 
in  similar  systems.  Specifically,  the  2D  packing  densities  of  anionic 
chains  or  rods  in  multi-lamellar  polyelectrolyte-membrane  com¬ 
plexes  (DNA  [25,26,28],  a-helical  polypeptides  [29])  and  polyelec¬ 
trolytes  adsorbed  at  membrane-solution  interfaces  (DNA  [30], 


polystyrene  sulfonate  [31])  were  found  to  increase  with  the  cat¬ 
ionic  lipid  concentration. 

However,  the  dependence  of  the  inter-TMV  spacing  on  the 
lipid-charge  density  that  we  observed  is  much  weaker  than  that 
reported  for  those  other  2D  polyelectrolyte-lipid  systems.  Defining 
/to  be  the  ratio  between  the  multiplicative  factors  for  the  concom¬ 
itant  increases  in  1  /d  and  in  the  cationic  lipid  fraction,  we  obtain 

1.2/20  =  0.06  (raising  XTAP  from  0.048  to  1.0  leads  to  a  ~20% 
increase  in  1  /d  or  ~rjl(r)  for  the  2D  assemblies  of  TMVs.  By  contrast, 
the  corresponding  /  values  are  higher  by  one  order  of  magnitude 
for  other  charged  rods:  /=  0.5-0.8  for  multi-lamellar  DNA-lipid 
and  polypeptide-lipid  complexes  [26,29], /~  1  for  DNA  adsorbed 
on  cationic  monolayers  [30],  and  f~  0.5  for  monolayer-bound 
polystyrene  sulfonate  [31]. 

This  apparent  discrepancy  again  may  be  explained  by  the  over¬ 
abundance  of  the  cationic  lipids  in  the  monolayer.  As  discussed 
above,  over  most  of  the  probed  range  of  lipid  compositions,  the 
charge  density  on  the  lipid  monolayer  should  far  exceed  that  on 
the  TMV  surface  (<~0.1  e/nm2).  While  the  counter-ion  release 
mechanism  implies  higher  adsorption  for  weaker  particle  charge 
and  higher  monolayer  charge  density  [14],  adsorption  is  eventually 
limited  by  the  electrostatic  (and  ultimately  steric)  repulsion  be¬ 
tween  adsorbed  particles.  In  addition,  TMV’s  large  size  and  surface 
curvature  implies  that  the  total  area  of  the  TMV-monolayer  con¬ 
tacts  should  be  much  smaller  than  the  total  interfacial  area.  For 
these  reasons,  the  interfacial  assembly  of  TMVs  is  expected  to  re¬ 
lease  only  a  small  fraction  of  the  counter-ions  that  are  originally 
associated  with  the  cationic  lipid  monolayer.  This  may  explain 
why  the  density  of  adsorbed  TMVs  does  not  strongly  depend  on 
the  lipid-charge  density.  By  contrast,  with  typical  surface  charge 
densities  of  0.2-1. 5  e/nm2  [27,29],  the  anionic  rods  are  more  highly 
charged  in  the  multi-lamellar  polyelectrolyte-lipid  complexes, 
whose  stability  at  low  ionic  strength  relies  on  close  matching  of 
the  charge  densities  on  the  polyelectrolyte  and  lipid-membrane 
surfaces  [27,29]. 

3 A. 3.  Adsorption  kinetics  affects  structural  order 

The  AFM  images  in  Fig.  3  further  reveal  that  relatively  large  do¬ 
mains  (>~1  pm)  of  aligned  TMVs  are  formed  on  the  monolayer  at 
XpAP  =  0.16  (Fig.  3A).  By  contrast,  TMVs  on  the  pure  DOTAP  mono- 
layer  (XTAP  =  1.0)  form  only  small  randomly-oriented  “clusters”  of 
five  or  fewer  rods  (Fig.  3C).  The  domain  sizes  observed  at  XTAP  =  0.4 
(Fig.  3B)  are  intermediate  between  the  above  two  cases.  This 
apparent  decrease  in  domain  size  with  increasing  lipid-charge  den¬ 
sity  is  qualitatively  consistent  with  the  GISAXS  results,  in  which 
higher  lipid-charge  density  results  in  a  reduction  in  the  modulation 


Fig.  3.  GISAXS  patterns  and  AFM  images  from  three  samples  under  Ca2+-free  10  mM  MES  buffers  with  three  different  lipid  monolayer  compositions:  XTAp  =  0.16  (A),  0.4 
(B)  and  1.0(C). 
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amplitude  of  S{qa)  (Fig.  SI),  indicating  a  greater  degree  of  posi¬ 
tional  disorder  in  the  TMV  assembly. 

The  dependence  of  structural  order  in  the  TMV  assemblies  on  li¬ 
pid-charge  density  may  be  the  consequence  of  the  adsorption 
kinetics  in  the  process  of  TMV  self-assembly.  Simple  estimates 
based  on  the  Grahame  equation  and  the  Gouy-Chapman  theory 
[32]  (Fig.  S2)  show  that  for  the  10  mM  MES  solution  at  pH  6,  raising 
the  DOTAP  fraction  from  5%  to  100%  is  expected  to  increase  the 
interfacial  potential  by  a  factor  of  ~3.5  [i j/(z  =  0)  ~  60  210  mV  > 

kBTle  =  26  mV],  leading  to  a  roughly  two-fold  enhancement  in  the 
potential  at  one  Debye  length  away  from  the  interface 
[i/^(z  =  4.8  nm)  ~  20  -►  38  mV].  Thus,  an  increase  in  the  lipid- 
charge  density  should  raise  the  adsorption  rate  by  strengthening 
the  electrostatic  lipid-TMV  attraction  near  the  interface.  The  faster 
TMVs  adsorbed,  the  less  time  would  be  allowed  for  the  TMVs  to 
freely  diffuse  and  reorient  themselves  after  adsorption.  This  in  turn 
should  induce  the  clustering  of  aligned  rods  at  more  sites  and 
cause  the  resulting  domains  to  be  smaller.  Therefore,  the  observed 
trend  in  the  domain  size  is  consistent  with  an  increase  in  the 
adsorption  rate  at  high  lipid-charge  densities. 

3  A  A.  Limitation  of  GISAXS  measurements  and  analysis 

While  we  did  follow  the  formation  of  TMV  assemblies  with  GI¬ 
SAXS  as  soon  as  TMVs  were  added  into  the  buffer  solution  (e.g., 
Fig.  IB),  our  data  did  not  conclusively  show  the  above  lipid-charge 
dependence  of  TMV  adsorption  kinetics.  In  Fig.  4A,  the  time  depen¬ 
dence  of  the  best  fit  r\\o  is  plotted  for  a  few  different  lipid  compo¬ 
sitions.  The  results  indicate  that  in  the  absence  of  Ca2+,  aligned 
domains  of  adsorbed  TMVs  started  to  form  almost  immediately 
after  injection.  However,  the  observed  time  evolution  of  q/a  did 
not  show  a  definite  systematic  trend  with  varying  lipid  composi¬ 
tion.  This  apparent  discrepancy  from  the  AFM  results  is  likely 
due  to  the  limited  sensitivity  of  GISAXS  to  adsorption.  The  above 
AFM-based  interpretation  about  the  lipid-composition  effect  on 
the  TMV  adsorption  rate  refers  to  the  earliest  stage  of  the  interfa¬ 
cial  assembly.  On  the  other  hand,  as  shown  in  Fig.  4A,  by  the  time 
the  first  set  of  analyzable  GISAXS  data  was  obtained,  the  local 
packing  density  had  already  reached  about  a  half  of  the  fi¬ 

nal  value.  In  view  of  this,  a  worthwhile  extension  of  the  present 
study  would  be  to  address  the  lipid-charge  dependence  of  the 
adsorption  rate  more  directly.  Doing  so  will  require  experimental 
methods  that  are  capable  of  finer  time  resolution  and  higher  sensi¬ 
tivity  to  adsorption,  such  as  quartz  crystal  microbalance  with  dis- 


t  (minutes) 

Fig.  4.  Time  dependence  of  the  best-fit  rj/a,  normalized  to  the  final  value  (q/c tf,  for 
three  different  lipid  compositions  under  10  mM  buffer  solutions  that  contain  (A) 
0  mM  and  (B)  100  mM  CaCl2. 


sipation  (QCM-D)  that  has  been  utilized  to  quantify  the  formation 
of  lipid  bilayer  on  a  substrate  [33]. 

With  the  decreased  domain  sizes  at  high  XTAP,  the  hard  rod 
model  that  approximates  the  TMV  assemblies  as  ID  hard-wall  fluid 
may  no  longer  be  adequate.  This  is  mainly  because  the  positional 
correlation  between  TMVs  within  each  domain  is  now  coupled 
with  their  orientations,  which  in  turn  are  influenced  by  the  pres¬ 
ence  of  TMVs  outside  of  the  given  domain.  Furthermore,  TMVs  that 
do  not  belong  to  any  ordered  domains  are  expected  to  scatter  X- 
rays  as  free  particles  (structure  factor  5  ~  1 ).  With  this  contribution 
in  the  observed  scattering  intensity,  the  contrast  of  the  TMV-TMV 
correlation  peak  would  be  reduced,  causing  the  data  to  deviate 
from  the  line  shape  expected  for  a  hard-wall  fluid.  The  inadequacy 
of  the  hard  rod  model  may  partly  explain  the  large  error  bars  on 
the  fitting  parameters  and  why  the  extracted  o  at  the  highest  DO¬ 
TAP  concentrations  falls  slightly  below  the  TMV  diameter  of  18  nm 
(Fig.  3A).  However,  this  complication  does  not  alter  the  main  con¬ 
clusions  about  the  qualitative  behavior  of  the  lipid-charge  depen¬ 
dence  of  the  interaction  between  TMVs,  namely,  high  XTAP  results 
in  smaller  average  spacing  between  TMVs  (higher  packing  density) 
but  poorer  positional  order.  These  conclusions  can  be  inferred  di¬ 
rectly  from  the  position  and  width  of  the  TMV-TMV  correlation 
peak  in  the  scattering  data  (see  Fig.  2A  and  Fig.  S3  in  Supporting 
information). 

3.2.  Effect  of  lipid  monolayer  charge  density  on  TMV  assembly  in  the 
presence  of  1 00  mM  Ca2+ 

Our  previous  study  for  XTAP  =  0.16  showed  that  the  presence  of 
Ca2+  ions  enhanced  the  structural  order  within  the  TMV  assemblies 
[10].  To  find  out  how  the  Ca2+  ions  would  influence  the  lipid- 
charge  dependence  observed  in  Section  1,  we  carried  out  a  similar 
set  of  experiments  at  various  values  of  XTAP  using  the  solution  that 
contained  100  mM  Ca2+  ions  (100  mM  CaCl2,  lOmM  MES, 
pH  ~  6.5).  The  observed  scattering  data  are  summarized  in 
Fig.  2B.  The  best  fits  based  on  the  sticky  rod  model  (solid  lines) 
reproduce  the  measured  scattering  profiles  reasonably  well.  The 
corresponding S(qxa )  curves  are  shown  in  Fig.  SI.  The  lipid-compo¬ 
sition  dependence  of  the  extracted  o  and  q\o  is  plotted  as  red  solid 
symbols  in  Fig.  2C. 

3.2 A.  The  screening  of  electrostatic  interactions  by  Ca2+  ions  reduces 
the  lipid-charge  dependence  and  growth  rate  of  the  TMV  assembly 

The  data  in  Fig.  2B  and  C  show  that  in  the  presence  of  Ca2+  both 
o  and  r\\o  are  essentially  independent  of  the  lipid  composition. 
Both  the  insensitivity  of  the  assembled  structure  to  the  cationic  li¬ 
pid  concentration  and  the  near  agreement  between  the  best-fit 
hard-wall  radius  o  (~19  nm)  and  the  TMV  diameter  (18  nm)  are 
consistent  with  strong  screening  of  the  electrostatic  inter-TMV 
repulsion  that  is  expected  for  the  100  mM  CaCl2  solution  (Debye 
length:  k~a  =  0.6  nm).  Note  that  in  contrast  to  the  Ca2+-free  case 
(Section  1),  ordered  TMV  structures  formed  on  the  pure  DOPC 
monolayer  in  the  presence  of  100  mM  Ca2+  ions.  This  may  be  due 
to  charging  of  the  monolayer  as  a  result  of  Ca2+  ions  binding  to 
the  phosphate  moiety  of  the  PC  headgroup.  Based  on  the  reported 
binding  constant  of  I<  =  12-20  M-1  [34],  we  estimate  the  surface 
charge  density  of  the  pure  DOPC  monolayer  to  be  as  high  as 
~20%  of  the  value  for  the  pure  DOTAP  monolayer  (Fig.  S4  in 
Supporting  information). 

Strong  screening  provided  by  Ca2+  ions  also  affects  the  growth 
of  ordered  domains.  As  shown  in  Fig.  4B,  the  time  dependence  of 
the  best-fit  q\o  indicates  that  the  assembly  develops  rather  slowly 
in  the  presence  of  Ca2+,  where  it  takes  2-3  h  for  the  local  packing 
density  {~rjl(j)  to  reach  ~50%  of  the  final  values.  This  should  be 
contrasted  from  the  Ca2+-free  case  (Fig.  4A),  for  which  r\\a  exceeds 
the  50%  mark  in  less  than  ~1  h.  The  presence  of  100  mM  Ca2+  is 
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expected  to  significantly  reduce  the  surface  potential  at  the  mono- 
layer-solution  interface  (by  a  factor  of  ~2;  see  Figs.  S2  and  S4)  and 
to  cause  the  electrostatic  TMV-lipid  attraction  to  be  much  shorter 
ranged  (kt1  =  0.6  nm  at  100  mM  Ca2+  vs.  4.8  nm  at  0  mM  Ca2+).  All 
these  effects  should  slow  down  the  TMV  adsorption  and  the 
growth  of  ordered  domains  and  contribute,  together  with  the 
Ca2+-induced  inter-TMV  attraction  to  be  discussed  next,  to  reduc¬ 
ing  the  influence  of  lipid  compositions  on  the  2D  assembly. 

3.2.2.  Attraction  between  TMVs  is  essential  in  forming  the  ordered 
assemblies  induced  by  Ca2+ 

As  noted  previously  [10],  a  key  result  of  the  above  analysis  is 
that  S{qx)  for  the  Ca2+-induced  assemblies  exhibit  relatively  sharp 
peaks  whose  magnitudes  decay  very  slowly  with  qx  (Fig.  SI).  This 
slow  decay  of  S(qx)  is  characteristic  of  the  sticky  rod  model  and 
is  found  to  be  essential  for  reproducing  the  data.  The  good  agree¬ 
ment  between  the  data  and  the  sticky  fluid  model  suggests  that 
the  overall  inter-particle  potential  for  lipid-bound  TMVs  includes 
a  short-range  attraction.  We  previously  speculated  that  this  in¬ 
ter-TMV  attraction  was  counter-ion-mediated  in  origin  [10].  This 
interpretation  is  indirectly  supported  in  the  present  study  by  the 
observed  insensitivity  of  the  Ca2+-induced  TMV  assemblies  to  the 
cationic  lipid  concentration.  It  should  be  noted  that  similar  obser¬ 
vations  were  previously  reported  for  multi-lamellar  DNA-cationic 
lipid  complexes,  for  which  high  concentrations  of  divalent  cations 
were  found  to  induce  close  packing  in  2D  arrays  of  parallel  DNA 
rods  [35].  In  that  study,  the  spacing  between  condensed  DNAs 
was  also  found  to  be  independent  of  the  cationic  lipid  concentra¬ 
tion,  and  this  behavior  was  interpreted  as  evidence  for  an  inter- 
DNA  attraction  that  was  mediated  by  correlations  between  the 
DNA-bound  divalent  cations.  Therefore,  the  similarities  between 
the  two  systems  suggest  that  an  analogous  Ca2+-mediated  attrac¬ 
tion  may  exist  in  the  2D  arrays  of  close-packed  TMVs.  The  origin 
of  the  inter-TMV  attraction  is  discussed  further  in  Section  3  below. 

3.3.  Ca2+  concentration  effect  and  nature  of  Ca2+ -induced  attraction 
between  TMVs 

In  order  to  gain  further  insights  into  the  effects  of  Ca2+  on  inter- 
TMV  interactions,  we  probed  the  structure  of  TMV  assemblies  at 


intermediate  Ca2+  concentrations  of  25  mM  and  50  mM.  At  Ca2+  con¬ 
centrations  above  150  mM,  aggregates  of  TMVs  formed  in  the  bulk 
solution  (see  Fig.  S5  in  Supporting  information).  The  in-plane  scat¬ 
tering  profiles  obtained  for  0, 25, 50,  and  1 00  mM  Ca2+  are  compared 
in  Fig.  5  for  two  lipid  compositions.  In  both  cases,  the  in-plane  corre¬ 
lation  peaks  shifted  to  higher  qr  values  as  the  Ca2+  concentration  was 
increased,  suggesting  a  decrease  of  the  average  inter-rod  spacing. 

These  results  reinforce  the  conclusion  from  Section  2  that  the 
Ca2+  ions  screen  the  electrostatic  repulsion  and  induce  attraction 
between  TMVs.  Fig.  5C  shows  that  the  presence  of  Ca2+  ions  signif¬ 
icantly  reduces  the  range  of  repulsion  between  TMVs.  As  the  Ca2+ 
concentration  is  increased  from  25  mM  (k:-1  =  1.2  nm)  to  100  mM 
(kt1  =  0.6  nm),  the  best-fit  value  of  o  decreases  from  ~20.6  nm 
to  ~19.0nm,  for  both  XTAP  =  0.048  and  0.3.  This  agreement  be¬ 
tween  the  two  values  of  XTAP  suggests  that  the  Coulomb  repulsion 
between  TMVs  is  screened  mainly  by  Ca2+  ions  and  not  by  the  cat¬ 
ionic  lipids  when  the  Ca2+  concentration  is  higher  than  25  mM.  On 
the  other  hand,  as  the  Ca2+  concentration  is  raised  from  0  to 
100  mM,  rjl<r  increases  by  a  factor  of  ~2,  more  or  less  monotoni- 
cally.  The  difference  between  the  qla  values  for  the  two  lipid  com¬ 
positions  again  vanishes  at  high  Ca2+  concentrations  (^50  mM), 
suggesting  dominant  Ca2+-induced  attraction.  In  order  to  distin¬ 
guish  whether  the  primary  role  of  Ca2+  in  inducing  close-packed 
structures  is  to  mediate  attraction  or  to  screen  repulsion,  and  to 
understand  the  origin  of  the  Ca2+-induced  attraction  between 
TMVs,  we  performed  the  following  additional  experiments. 

3.3.1.  Screening  alone  cannot  induce  ordered  TMV  assemblies 

We  prepared  TMV  samples  following  the  procedures  described  in 
Section  1  for  three  lipid  compositions:  XTAP  =  0.6,  0.8,  and  1.0.  After 
allowing  TMV  assemblies  to  develop  for  -42  h  and  GISAXS  patterns 
(Fig.  6A)  confirmed  the  formation  of  TMV  assemblies,  we  rinsed  the 
sample  cell  with  the  Ca2+-free  buffer  several  times  to  remove  the  ex¬ 
cess  TMVs  in  the  solution.  We  then  added  a  small  amount  of  concen¬ 
trated  Ca2+  solution  to  the  remaining  solution  to  achieve  the  final 
Ca2+  concentration  of  -400  mM.  After  another  12  h  of  incubation, 
GISAXS  patterns  were  again  recorded  (Fig.  6A  and  B).  The  corre¬ 
sponding  best-fit  values  of  rj/a  are  plotted  in  Fig.  6C. 

For  all  three  lipid  compositions  probed,  the  addition  of  Ca2+  re¬ 
sulted  in  substantial  changes  in  the  GISAXS  patterns.  The  in-plane 


Fig.  5.  In-plane  scattering  profiles  from  TMV  assemblies  adsorbed  to  lipid  monolayers  with  XTap  =  0.3  (A)  and  0.048  (B),  respectively,  under  MES  buffers  containing  100  mM 
(squares),  50  mM  (circles),  25  mM  (triangles)  and  0  mM  (stars)  Ca2+.  The  data  are  shifted  vertically  for  clarity.  The  solid  lines  are  the  best-fits  based  on  the  hard  rod  model  for 
the  Ca2+-free  case  and  the  sticky  rod  model  for  the  non-zero  Ca2+  concentrations.  The  corresponding  best-fit  values  for  the  radius  of  hard-wall  potential  a  and  the  packing 
density  of  rods  rj/c t  are  shown  in  (C). 


S.  Wang  et  al. /Journal  of  Colloid  and  Interface  Science  358  (201 1 )  497-505 


503 


Fig.  6.  The  in-plane  scattering  profiles  from  three  samples  withXTAP  =  0.6  (squares),  0.8  (circles),  and  1.0  (triangles)  (A)  on  the  lipid  monolayers  under  a  Ca2+-free  MES  buffer, 
and  (B)  from  the  same  samples  at  12  h  after  increasing  Ca2+  concentration  to  100  mM.  The  insets  are  GISAXS  patterns  from  the  sample  with  the  XTAP  =  0.6  lipid  monolayer. 
Solid  lines  are  the  best  fits,  with  the  best-fit  packing  density  r\\o  plotted  in  (C).  Note  that  the  best-fit  parameters  for  the  data  collected  after  increasing  the  Ca2+  concentration 
are  different  from  those  presented  in  Fig.  2C  for  samples  with  the  same  final  Ca2+  concentration.  This  is  most  likely  because  in  the  present  case  the  overall  density  of  adsorbed 
TMVs  is  lower,  without  free  TMVs  to  be  adsorbed  after  the  increase  in  Ca2+  concentration. 


correlation  peaks  became  more  prominent  and  sharper,  indicating 
that  the  TMV  arrays  became  better  ordered.  The  qr  positions  of  the 
first-  and  second-order  peaks  (0.032  A-1  and  0.064  A-1)  that  devel¬ 
oped  after  the  addition  of  Ca2+  are  almost  independent  of  the  lipid 
compositions  (Fig.  6B)  and  correspond  to  an  average  spacing  of 
~20  nm  between  TMVs.  Given  the  AFM-derived  spacing  of  32- 
35  nm  before  adding  Ca2+  (Section  1),  the  Ca2+-induced  reduction 
in  the  inter-TMV  spacing  was  12-15  nm  (~40%).  These  results 
demonstrate  that  the  Ca2+-induced  close  packing  of  lipid-bound 
TMVs  is  brought  about  by  an  attractive  inter-TMV  interaction. 
Had  the  added  Ca2+  ions  only  screened  the  electrostatic  repulsion 
between  TMVs  without  introducing  attraction,  the  packing  density 
in  the  TMV  assemblies  (or  average  spacing  between  TMVs)  would 
have  remained  the  same  because  the  solution  contained  no  excess 
TMVs  to  adsorb  to  the  lipid  monolayer. 


3.3.2.  Attraction  between  TMVs  has  a  contribution  that  depends  on 
counter-ion  valence 

As  pointed  out  earlier  (Section  2),  the  success  of  the  sticky  rod 
model  in  describing  the  Ca2+-induced  structure  strongly  suggests 
that  short-range  attractive  forces  must  exist  between  the  lipid- 
bound  TMVs  to  facilitate  the  formation  of  close-packed  assemblies. 
The  most  probable  inter-rod  spacing  is  determined  by  the  balance 
between  the  repulsive  and  attractive  inter-particle  forces.  Once  the 
repulsive  contribution  is  reduced  by  more  effective  screening  of 
the  surface  charges  on  TMVs,  the  attractive  contribution  is  ex¬ 
pected  to  draw  the  TMV  rods  closer  together  in  the  assemblies. 
This  attractive  force  between  like-charged  TMVs  should  partly 
arise  from  the  van  der  Waals  interaction,  which  is  thought  to  be 
responsible  for  the  attraction  between  colloid  particles  in  solution 
[36,37]  as  well  as  the  formation  of  TMV  bundles  in  bulk  solutions 
at  high  ionic  strength  [15].  As  previously  discussed  in  Ref.  [10],  part 
of  the  attraction  could  also  be  mediated  by  counter-ions,  which  is 
also  believed  to  be  responsible  for  the  lateral  aggregation  of 
charged  rods,  such  as  DNAs  [35]  and  other  filamentous  viruses 
[38],  under  high  concentrations  of  divalent  cations.  While  coun- 
ter-ion-mediated  attraction  is  expected  to  depend  on  the  valence 
of  the  counter-ion,  the  van  der  Waals  force  is  not.  We  therefore 


examined  the  effect  of  the  cation  valence  on  the  structure  of 
TMV  assemblies. 

We  prepared  TMV  samples  for  two  lipid  compositions, 
XTAp  =  0.3  and  0.048,  under  pH  6.0,  10  mM  MES  solutions  contain¬ 
ing  300  mM  NaCl,  which  have  the  same  ionic  strength  as  the 
100  mM  CaCl2  solution.  The  final  scattering  patterns  (Fig.  7) 
showed  clear  scattering  peaks  indicating  in-plane  ordering.  The 
in-plane  scattering  data  are  well  described  by  the  sticky  rod  model, 
suggesting  that  attractive  forces  also  exist  between  TMVs  in  the 
presence  of  Na+  ions.  However,  the  extracted  hard-wall  radius  of 
~22  nm  is  significantly  larger  than  the  value  of  ~19  nm  obtained 
for  Ca2+-induced  structures  (Section  2).  Moreover,  the  extracted 
stickiness  parameter  (see  Experimental  and  Ref.  [10])  of  a  ~  0.3 
for  Na+  is  significantly  lower  than  the  average  value  of  a  ~  1  for 
Ca2+.  Therefore,  the  van  der  Waals  force,  which  is  expected  to  be 
independent  of  the  type  of  cations  in  the  solution,  cannot  solely  ac¬ 
count  for  the  entire  attraction  between  TMVs  in  the  presence  of 
Ca2+  ions. 

On  the  other  hand,  the  counter-ion-mediated  attraction  be¬ 
tween  charged  rods  is  thought  to  arise  from  density  fluctuation 
in  the  counter-ions  that  are  condensed  on  the  rods  [23].  The  con¬ 
densation  of  counter-ions  is  expected  to  occur  when  the  linear 
charge  density  on  the  rod  is  sufficiently  high  that  the  electrostatic 
potential  energy  overcomes  entropy  [23].  The  rod’s  linear  charge 
density  is  often  measured  by  the  unitless  parameter  £  =  lblb ,  where 
lb  is  the  Bjerrum  length  and  b  is  the  length  of  the  rod  per  unit 
charge.  Counter-ion  condensation  occurs  when  £  exceeds  a  critical 
value  £c.  The  resulting  condensation  then  normalizes  the  nominal 
linear  charge  density  so  that  £  =  £c.  For  ideal  cylinders  [39],  £c  scales 
with  1/z,  where  z  is  the  valence  of  the  counter-ions.  Therefore, 
counter-ion  condensation  is  more  likely  to  occur  for  ions  of  higher 
valence,  resulting  in  a  lower  renormalized  linear  charge  density  on 
the  rod,  and  consequently  weaker  electrostatic  inter-rod  repulsion 
and  stronger  attraction  that  is  mediated  by  the  density  fluctuations 
in  the  condensed  counter-ions.  This  expected  trend  is  consistent 
with  our  observation  that  Ca2+  ions  introduce  stronger  attraction 
between  TMVs  than  Na+  ions  do,  suggesting  that  counter-ion-in- 
duced  attraction  does  play  a  role  in  the  formation  of  TMV  assem¬ 
blies  in  the  presence  of  Ca2+  ions. 
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Fig.  7.  In-plane  scattering  profiles  from  the  final  GISAXS  pattern  (inset)  of  two  samples  with  (A)Xtap  =  0.3  and  (B)XTAP  =  0.05  under  pH  6.0, 10  mM  MES  solutions  containing 
300  mM  NaCl.  The  open  circles  are  the  data.  The  solid  lines  are  the  best-fits  using  the  sticky  rod  model. 


4.  Conclusion 

We  have  studied  the  structure  of  2D  assemblies  of  TMVs  ad¬ 
sorbed  on  lipid  monolayers  of  DOPC/DOTAP  mixtures  at  various 
mole  fractions  of  the  cationic  DOTAP  and  under  buffer  solutions 
containing  different  types  and  concentrations  of  cations.  These  fac¬ 
tors  influence  the  electrostatic  interaction  between  TMVs  as  well 
as  that  between  TMV  and  the  lipid  monolayer,  which  in  turn  im¬ 
pact  both  the  behavior  of  TMV  adsorption  onto  the  lipid  monolayer 
and  the  final  structure  of  the  TMV  assemblies. 

In  general,  the  2D  assembly  of  TMVs  is  found  to  be  slower  under 
solutions  of  higher  ionic  strengths,  which  can  be  attributed  to 
screening  of  the  electrostatic  attraction  between  the  TMVs  and 
the  lipid  monolayer.  At  low  ionic  strength  (10  mM  MES  only),  the 
reduced  screening  causes  the  assembly  to  occur  more  quickly.  Un¬ 
der  these  low-salt  conditions,  increasing  the  fractions  of  cationic 
lipids  in  the  monolayer  resulted  in  poorer  structural  order.  This 
observation  can  be  attributed  to  accelerated  adsorption  of  TMVs 
at  high  lipid-charge  densities,  which  we  interpret  as  the  conse¬ 
quence  of  higher  surface  potential  at  the  lipid  monolayer.  Besides 
the  disordering  effect,  higher  DOTAP  concentration  in  the  lipid 
monolayer  also  resulted  in  higher  TMV  packing  density,  possibly 
due  to  the  partial  screening  of  electrostatic  repulsion  between  ad¬ 
sorbed  TMVs  by  mobile  counter-ions  within  the  lipid  monolayer 
(i.e.  DOTAP).  Modification  of  the  interaction  between  adsorbed 
TMVs  by  the  ions  in  the  bulk  solution  was  found  to  dramatically  af¬ 
fect  the  TMV  assemblies,  yielding  denser  and  better-ordered  struc¬ 
tures  with  increasing  concentration  of  Ca2+  ions.  The  improved 
structural  order  may  be  attributed  to  increased  screening  of  the 
Coulomb  repulsion  as  well  as  counter-ion-induced  attraction  be¬ 
tween  TMVs. 

Our  results  revealed  the  complexity  of  the  electrostatic  interac¬ 
tions  that  are  involved  in  the  2D  assemblies  of  TMVs  on  oppositely 
charged  lipid  monolayers.  While  we  were  able  to  resort  to  simpli¬ 
fied  models  to  quantify  the  final  structure  in  the  TMV  assemblies, 
the  initial  behavior  of  TMV  adsorption  cannot  be  inferred  from  the 
X-ray  scattering  data  because  ordered  assemblies  have  already 
formed  by  the  time  the  scattering  intensity  is  sufficiently  high  to 
be  analyzed.  The  ID  liquid  model  used  in  our  data  analysis  should 
also  break  down  when  the  TMVs  have  not  yet  formed  large,  or¬ 
dered  domains.  Our  data  cannot  answer  definitively  whether  the 
attractive  interaction  between  TMVs  in  the  presence  of  divalent 


cations  is  due  to  the  condensation  of  these  ions  onto  the  TMVs. 
Nevertheless,  our  results  provide  strong  evidence  that  the  inter- 
TMV  attraction  has  a  component  that  depends  on  the  counter¬ 
ion  valence.  The  role  of  counter-ion  condensation  may  be  further 
clarified  by  X-ray  scattering  measurements  from  TMVs  in  bulk 
solution  in  the  presence  of  high  Z  counter-ions. 
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